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ABSTRACT 


The Naval Postgraduate School simulation model, QDYN&6, was used to 
examine sputtering of nitrogen from the (100) faces of single crystals 
of molybdenum and tungsten. The nitrogen placement was varied, = and 
analyses were conducted on the sputtering cross sections of the 
nitrogen. The cases where the adatom was directly hit by the incident 
1on, or sf it was sputtered due to the collision cascade process, were 
analyzed separately. The simulations were conducted to compare the 
results with Winters’ recent work, and to build upon the efforts of ear- 
lier studies compieted at the Naval Postgraduate School. lt was found 
that placement of nitrogen at 0.245 A from the surface of molybdenum 
resulted in cross sections similar to those found by Winters. The 
effect of the mass of the substrate was verified, in that a substrate of 
greater mass results in a higher sputtering cross section. This agreed 
with Winters’ findings, and conflicted with earlier conclusions of past 
theses. The adatoms apparently reduce the momentum available to create 
collision cascades, reducing the sputter yield ratio of the substrate 


when the ions directly hit the adatoms. 
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F BACKGROUND 


A, HISTORICAL OVERVIEW 

When a surface is bombarded by a beam of energetic particles under 
proper conditions, surface damage effects can be observed. This damage 
can mamfest itself a number of different ways; as surface damage in the 
form of opits, blisters, and cones, and in the ejection of atoms from the 
target. This eyection of atoms during bombardment is called sputtering. 

Sputtering was first discovered in 1853 by Grove [Ref. 41} when he 
observed the disintegration of the cathodes in glow-discnarge tubes. He 
noted that the cathode material was deposited on the glass walis of the 
tubes. He called this process cathode sputtering. This erosion of the 
cathode was also noted by Plucker [Refs. 2-4]. Gassiot [Ref. 5] and 
Faraday also reported similar observations. 

Fifty years passed before any explanation for these phenomena was 
proposed. In 1902 Goldstein [Ref. 6] presented evidence that the sput- 
tering effects reported by Grove and Faraday were caused by positive 
atoms of the discharge impacting on the metal cathode of the tubes. 
Seven years later, in 1909, Stark proposed two models in an attempt to 
explain sputtering. The "hot spot model" considered the sputtered atoms 
to be the result of evaporation of target material from a small surface 
region due to localized heating by the i10n beam. The "collision model” 
proposed that sputtering events were the result of a series of binary 
collisions imtiated by a single ion {[Ref. 7). These models served to 
explain the experimental results at the time. 
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Later, tw 14921, Thompson mike. 8} contributed to the field by 
suggesting that the atomic eyection was caused by the release of radia- 
tion as the ton struck the target. in the following years, a consider- 
able Senne of theoretical and experimental work was done. Bush and 
Smith (Ref. 9] in 1922 attempted to describe the eyection of atoms as 
the results of the expansion of gas adsorbed by the target material. The 
following year, Kingdon and Langmuir [Refs. 10, ‘i1) conducted an exper- 
iment that suggested a momentum transfer ejection mechanism for sputter- 
ing. They bombarded thoriated tungsten with tons in a glow discharge 
tube. This was a special case of sputtering, where the thin surface film 
of thorium sputtered rather than the tungsten substrate. 

In 1926, von Hipple and Blechschmidt§ [Refs. 12-15] proposed a 
theory that described sputtering 35 an evaporation of the surface atoms. 
Through spectroscopic techniques, von Hipple found that some of the 
sputtered atoms were tn an excited state. He made further refinements 
on Stark's hot spot model, and made the first attempt to formulate a 
Sputtering theory on the basis of local heating. 

Lamar and Compton [Ref. 16}, in 1934, publisned “A Special Theory 
of Cathode Sputtering," which led to the "thermal spike" concept. They 
suggested that binary collision processes were dominant «tn light-ion 
sputtering, where local evaporation predominates for heavy-ion sputter- 
Ing. The thermal spike concept was based on momentum transfer between 
the bombarding ion and the lattice atoms. The theory suggested that a 
long hved high temperature volume persisted in the target after the 


collision cascade was completed. 


One of the major problems facing experimentalistSs was the iack of 
ability to reproduce results. Penning and Moubis in 1940 publisned = the 
results of their studies of the effect of pressure on sputtering yield 
[Ref. 17], (the sputtering yield is defined as the ratio of "sputtered" 
target atoms per incident ion). They found that the sputtering yield 
was reduced with an increase in- pressure. Collisions between the 
eyected atoms and the surrounding gases at the higher pressures resulted 
in the backscattering of eyected atoms back onto the target surface, 
reducing the sputtering yield. When the background pressure was Kept 
below 107° Torr, they were able to obtain reproducible results for ion 
energies in excess of 500 eV. 

Keywell made an attempt to formulate Stark’s collision model in 
terms of a neutron transport model originally developed for nuciear work 
(Ref. 18). Wehner [Refs. 19 - 21], one of the major contributors to the 
modern understanding of sputtering, began publishing his” findings itn 
1954, He discounted the evaporation model, and presented strong é€vi- 
dence for a momentum transfer process. He demonstrated the effects of 
crystal structure on the yield, and it became apparent that local heat- 
ing alone couid not account for the effects of sputtering. His observa- 
tions revived interest in collision theory. One of his mayor contribu- 
tions was the discovery of the anisotropic ejection patterns from mono- 
crystalline targets, the now well Known “Wehner spot pattern." 

In 1956 Harrison [Ref. 22] applied statistical methods to sputter- 
ING. He developed a theory involving the interaction of two distribu- 


tion functions; one for the crystal lattice and one for the ion beam. 


He utilized transport theory and introduced the 1Gea of Gea Colas 
cross sections to explain sputtering from amorphous materials. 
Sputtering, to this point, had been regarded as an annoying mani- 
festation, one that eroded cathodes and filaments, contaminated plasmas, 
and was a general nuisance. There was great difficulty in obtaining 
reproducible experimental results. The theories proposed had not com- 
pletely described the phenomena, and indeed have not yet today. Before 
continuing with the discussion of the current trends in sputveraing 
theory and experimentation, 2a brief discussion of the currently 


"accepted" concepts of sputtering is in order. 


B, PHYSICAL UNDERSTANDING 

Sputtering is the eyection of atoms from a target when bombarded by 
energetic ion projectiles. The incoming ion collides with atoms in the 
bulk of the target material, transferring part of its Kinetic energy 
and momentum to the target atoms. If the energy transferred to an atom 
i$ greater that the binding energy at the lattice site, then 42 primary 
recoil atom is created. This recoil atom may then collide with otner 
atoms in the lattice, creating a collision cascade, transferring energy 
throughout the material. A surface atom mey be ejected (sputtered) if 
the normal component of the energy transferred to it is greater than the 
surface binding energy of the lattice [Ref. 23]. Figure 1 illustrates 
three concepts of the sputtering process. Figure 1 (a) shows the inter- 
action of the incident ion beam with the surface layer of the tanger: 
Figure 1 (Dd) illustrates the concept of the “thermal spike," where a 


high ltocal temperature results in the evaporation of surface atoms. 
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Finally, Figure 1 (c) ijlustrates the collision cascade PROCESS, Which 


is the current conceptual model for sputtering [Ref. 24), 


Ion Beam 


+> 7 val 
Y 


(a) (b) (sep) 


(a) Sputtering of surface atoms aS a result of the inter- 
action of the incident ion beam and the target sur- 
face. 


(Db) Surface atoms evaporated as a result of a thermal] 
spike, 


(C) Sputtering from a_i collision cascade, where energy is 
directed toward the surface from multiple binary 
COMiSions, [Ref. 24) 


Figure 4, Models of the sputtering Process. 


The mechanisms for the sputtering of multicomponent materials 
(targets consisting of more than one atomic component) are more compli- 
cated. One of the simpler cases 1S the situation where a thin layer of 
a third component is placed on the surface of the thick target material. 
This iS simiar to the case examined by kingdon and Langmuir in Refer- 
ence 11, In this case, the added atoms on the surface are called adatoms 
and the bulk of the material ts called the substrate [Ref. 25}. Examples 
of this situation were covered by Garrison, Winograd, and _ Harrison 
fRef. 26} for oxygen on copper, and by Winters and S:gmund for nitrogen 
on tungsten [Ref. 27}. 

Winters proposed three mechanisms for the sputtering of the surface 
atoms, illustrated below in Figure e. Layer (3) is the thin Jayer of 
adatoms on the surface of the thick target (2). An ion (1) bombards the 
target surface, In Figure 2(a), the ion hits the adatom, and the adatom 
is reflected off the substrate, either directly, or after it penetrates 
the target slightly. Figure e(b) wWustrates the case where the ton 
does not hit the adatom directly, but penetrates the target, and is in 
turn reflected itself. The reflected ion then hits the adatom, which is 
then sputtered off. Figure e(c) shows the case where the ion causes 2n 
outward collision cascade, much like that illustrated in Figure i(c). 
The outward flux of sputtered substrate in turn sputters the adatom. 

The basic processes associated with sputtering are simiar to those 
causing radiation damage in the bulk of a= solid. Sputtering, however, 
usually involves atoms in the distorted surface layers (selvage) of the 
target material. The term, "selvage," is derived from the term “"“sel- 


vedge,” which iS a narrow band woven such that the edge will not 


2 


unravel. This analogy was extended to the distorted layers of the sub- 


Strate material at its surface (Ref. 28). 


(a) 





(1) 


. 7 


cep 





Figure 2. Three Proposed Mechanisms for the Sputtering of a Thin 
Layer on a Thick Target [Ref. 27]. 

Sputtering events take place in conditions far from thermal} equili- 
Drium, and are not evaporation of the material. several factors can 
influence the sputtering yield, including incident ion energy, angle of 
incidence, i0n type, target material, target crystal orientation, the 


presence Of adsorbed molecules on the surface of: the target, ambient 
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pressure, and numerous other’ factors. In order to obtamn reproducible 
experiments, the following conditions must be met {Ref. 29]: 
1. The target surface must be clean, that is free of contaminants in 
the form of adsorbed gases, lubricants from the vacuum pump, 


cleaning solvents, etc. 


e. The gas pressure must be such that the mean free path of 10ns and 
sputtered atoms is large. 


oy The ion current density must be mgh and the background pressure 
low so that formation of surface layers is prevented during the 
experiment. 

4, The ions must strike the target at a Known angle. 

= The energy spread of the incident beam must be_ small. 

sy The ionizing conditions in the ion source should be such 4s to 
minimize the production of multiply charged species; the ions 


must be unmformly charged and mass separated. 


Te The lattice orientation of monocrystalline targets must be Known. 


The field of sputtering is mch with terms used to describe the 
various observed and theoretical aspects of sputtering. Specific terms 
will be  iwtroduced as necessary to describe certain events, but two 
terms do merit mention here. Transmission sputtering is the ejection of 
atoms from the rear of a thin. target. This occurs when sufficient 
energy is transported throughout the target to allow atoms to overcome 
the potential energy binding it to the target, allowing it to escape. 
Back-sputtering is the more familiar form, where atoms are eyected from 
the surface of the target material. Two other terms widely used, espe- 
cially when discussing the sputtering of multicomponent materials are 
physical sputtering and chemical sputtering. Physical sputtering 
involves a transfer of Kinetic energy from the incident particle to the 


atoms in the target, and the subsequent ejection of the atoms. Back 
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Sputtering and transmission sputtering are manifestations of physical 
sputtering. CieiiegmeopuLtenring = = 15 tne result of a chemical reaction 
induced by the bombarding particle, which produces an unstable chemical 


compound on the surface of the target [{Ref. 30}. 


©, THEORETICAL AND EXPERIMENTAL DEVELOPMENTS (1960 170 PRESENT) 

In 1962, Wehner and others [Ref. 31] published data on the =  sput- 
tering yields of metals and semiconductors in the 100-600 eV range. 
Later in 1966, Wehner [Ref. 32], in a report for Litton Systems pub- 
lished detailed results covering years of low energy sputtering 
research. 

Silspee [Ref. 33], in an effort to account for the angular distri- 
bution of eyecta, the "Wehner spot patterns," proposed a focused colli- 
sion model, which allowed the transport of momentum in crystals along 
preferred directions. Experimental results had indicated that the 
Sputtering yield for. single crystal targets were dependent upon the 
crystallographic orentation of the target and the ineident ion beam. 
This seemed reasonable, considering the "holes" opserved in crystal 
models of when viewed from different planes. Focusing can be seen to 
contribute to the development of collision cascades within the target, 
it could not explain the Wehner spot patterns. 

Sigmund and Lehmann [Ref. 34] proposed an alternative model, based 
on Boltzman transport equations, requiring the target surface to have an 
ordered surface. In this model, an atom in the collision cascade 


would sputter if its Kinetic energy component normal to the surface was 


1S: 


sufficient to overcome the surface potential barrier. Thompson [Ref. 
35] proposed another model to account for the spot patterns, in which 
the surface attraction for the escaping atom causes a refraction of its 
velocity vector away from the normal, resulting in a distortion of the 
angular distribution of eyecta. The theories of Sigmund and Thompson 
predict the sputtering yields of amorphous or polycrystalline targets 
fairly well, but do not accurately reflect the experimental results from 
the ordered surfaces of single crystal targets. 

Numerous theories have been proposed, but none succeeds in fully 
describe all aspects of sputtering phenomena. Harrison [Ref. 36} in 2 
recent review, stated that there were currently no less than seven types 
of sputtering theories in the ‘literature. The statistical theories of 
Thompson [Ref. 35)) and Sigmund [Ref. 34} predict the yield relatively 
well, and provide information on the eyected atom energy § distribution 
function, and ejected atom angular distribution for  polycrystalime 
targets. Kelly (Ref. 37) classifies sputtering processes according to 
time scales. 

During this time, computer simulations were used to atin sputter- 
ing. The models and codes now in use have been evolving for the past 30 
years. The next section will present some basic philosophy of computer 


Simulations, with applications to sputtering research. 


D). COMPUTER SIMULATIONS 
In 1960 a new tool was added to the scientific arsenal, aimed at 
examining sputtering. Gibson, Goland, Milgram and Vineyard {[Ref. 38) 


published a report of the first computer simulation examining radiation 
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damage in a material. They simulated metallic copper and studied radia- 
tion damage events at low and moderate energies, up to 400 eV. In the 
model, one atom was given an arbitrary Kinetic energy and direction of 
motion, stmulating itsS having been struck by an energetic particle. 
This was one of the first published accounts of a computer simulation to 
study such -events. 

The Nhmgh operating speed of computers makes them the natural chorce 
for processing the numerous calculations required in a numerical = analy- 
sis. The use of a computer simulation frees one from the constrants of 
general theories, and allows one to examine the basic physics of the 
system, to see how theory and experiment interact. In the words of 
Harrison {Ref. Somape 4): 

A simulation is not a theory: it tS a mathematical tool which is used 
to test the fundamentals of a theory. The computer can model 4 
system with a minimum set of physical assumptions. This inherent 


simplicity helps to elucidate complex problems like sputtering: 
simulations develop ideas which can then be exploited by both 


experimentalists and theorists. 


Computer simulations fall into two general categories: time-step 
models and event-store models. An event-store program moves from event 
to event, skipping the tntervening§ time. It maintains a lst of poten- 
tral future events (hence the name "“event-store"), and checks and 
updates the lists to determine which event happens next. This model 
works well when the model is well understood, and the events are well 
separated in time. A time-step program carries the model forward for a 
Short period of time, computes everything that happens to the system tn 
that time interval, updates, and then continues on. This program 15 


most useful when several things happen simultaneously. The time-step 


ey 


program tends to be shorter tnan the event Stone “program, but (emai 
runs slower. 

The event-store model used in sputtering research is hased on the 
binary collision approximation (BC) [Ref. 39}. The assumption made is 
that each particle only interacts with one other particle at a time, and 
this other particle 1S usually assumed to be = stationary. These models 
are inherently linear calculations [Ref. 4Q}. 


The time-step model used in ths investigation is based on simulta- 


neous multiple interactions (Ml). Newtons laws, usually expressed in 
Hamiltonian form, are numerically solved for many particles. The Ml 


model used in this thesis, QDYN86, is the latest revision of the Mi 
program developed and used at the Naval Postgraduate School, and Penn- 
sylvania State University. QDYN86 is a full-lattice simulation, which 
models the dissipation of an incident ion’s momentum in a_e single crystal 
target, using classical mechanics. This program can generate different 
surfaces of several crystal structures. Adatoms can be placed on the 
surface of the crystal to simulate reacted surfaces. Specific details of 
the simulation will be presented m a later section of this” thesis. 

As can be seen, a considerable effort has been expended to study 
the mechanisms of sputtering. There 18 8 practical reason for the 


interest, beyond the lure of pure research. 


Ey APPLICATIONS 


Sputtering was long regarded as an undesired and lIttle-understood 


effect (Refs. 41, 42}. it destroyed cathodes and grids in gas discharge 


lke 


tubes, and contaminated plasmas and the surrounding walls (Ref. 43). 
There was great concern about damage to spacecraft and satellites from 
Sputtering. Efforts to understand and reduce these effects provided 
mayor impetus to the study of sputtering. 

However; sputtering, in the form of the controlled removal of 
Surface atoms from a target, is becoming especially important in its own 
right. The ability to precisely control an ion beam, and to remove 
atoms from surfaces 1S a very important tool in research, and in the 
manufacture of minature components. 

sputter ion sources can be used for cleaning of surfaces, by remov- 
ing adsorbed surface molecules to a degree that is impossible to achieve 
chemically or mechanically. The atoms sputtered from a surface can be 
analyzed in a mass spectrometer, giving information about the surface 
composition. One of the more important commercial applications of 
Sputtering if the deposition of thin films on a large variety of sub- 
Strates, especially useful in the manufacture of microelectronics. 
sputtering is used in micromachining and depth profiling of thin’ films 
fRef. 23). 

There are applications of erosionally modified surfaces in the 
fabrication of optical, magnetic and surface acoustical technologies. 
Areas include grating fabrication, magnetic bubble technology, ion 
polishing, and reactive ion-beam etching [Ref. 44}. Research is being 
conducted to examine possible biomedical applications of sputtering. 
These include surface modification of biomedical materials, using sput- 
tering for pathological discrimination, and with applications to 


implants and prostheses [Ref. 45}. 
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One of the first commercial applications of sputtering was the use 
of sputtering for the deposition of solid film lubricants (Ref. 46]. 
Sputtering had the advantage that it allowed for the deposition of 4 
variety of materials, on a large variety of substrates. The sputtered 
films are very dense, and strongly adherent. These properties are 
particularly useful for corrosion resistance and _ lubrication. 

The list of possible applhcations could <oitiiewe! The ability to 
understand and control the surface characteristics of a material to the 
atomic Jevel 18 a very powerful tool. The mayority of the work, espe- 
cially in regards to the applications, is the result of considerable 
experimental effort. if one could better understand the mechanisms 


underlying the observed manifestations, then experiments 


designed, and ‘applications could be more quickly 


this end that this research is conducted. 
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realized. 


could be better 


It is toward 





I. OBJECTIVES 


A, PREVIOUS SEPRORTS 

The interaction of nitrogen with the (100) plane of tungsten has 
been widely examined. It was considered to be a useful system on which 
to base the development of adsorption Kinetic and dynamic models 
{Ref. 47}. Winters [Ref. 48], in 1982, published a paper comparing the 
sputtering of chemisorbed nitrogen from polycrystalline targets of 
molybdenum and tungsten. His experiments showed that the mtrogen 
sputtering yield tended to increase as the atomic weight of the target 
substrate increased. 

This paper of Winters provided the basis for two Masters theses at 
the Naval Postgraduate School. In 1983, Meyerhoff [Ref. 49] used an 
earher version of Harrison’s computer simulation sputtering model to 
Study sputtering from nitrogen reacted (100) surfaces of tungsten and 
molybdenum targets. His g0al was tO compere the simulation results 
with the experimental results recently published by Winters in Refer- 
ence 48. Meyernoff concentrated on the reported mass effects. He 
concluded that while there may be a mass relation, the distance of the 
adatom from the substrate had a much more profound effect on the sput- 
tering yield. 

In 1986, Webb [Ref. 50} sought to further examine the adatom pla- 
cement problem. Webb again used (100) surfaces Of monocrystalline 


tungsten and molybdenum reacted with mtrogen and bombarded with argon 


YM 


ions. He narrowed the range of the placement of the adatoms, and nis 
results indicated that the adatoms were positioned above the (100) sur- 
face of the target plane. 

In all these cases, the nitrogen adatoms were assumed to be in the 
four-fold position, as shown below in Figure 3. Meyernoff assumed the 
nitrogen atom was slightly above the surface of the target plane, in a 
position of equal distance from its five nearest neighbors. Webb exa- 
mined the nitrogen sputtering yield with the nitrogen placed at three 


different heights for molybdenum and two heights for tungsten. 


—_ Adatom ireaeold 
~~ Position 
“ee, 






(100) SUmiaee 






Body- 
Centered t } 
Atom Se 


Figure 3. Nitrogen Adatom Placed in Four-fold Position in BCC 
Unit Cell. 

The basis for the four-fold placement assumption is well estab- 
ished in the literature. An earher work by Clavenna and schmidt 
examined the interaction Of No with W(100) [Ref. 5/4). They studied the 
binding states and condensation and desorption Kinetics of nitrogen on 


W(100) using flash desorption spectrometry. They felt that the 


Ze 


nitrogen was in the four-fold position, Out were uncertain exactly 
where the nitrogen was placed vertically in relation to the target sur- 
face. Adams and Germer [Refs. 52, 53] also exammed the adsorption of 
nitrogen a tungsten. They used a combination of Low Energy Electron 
Oiffraction (ERED); flash-desorption mass-spectrometry and contact 
potential measurements in their experiments. They refined thetr exper- 
iment by using the ION-isotope of nitrogen, in order to distinguish 
between nitrogen and CO present in the background gas. They felt that 
since the mitrogen atom was half the size of the tungsten atom, that it 
would sit in the well formed by the surface atoms of the substrate. 
Griffiths, Kendon, King and Pendry [Ref. 54] examined target atom 
displacement due to the introduction of the nitrogen adatoms. They 
performed as series of LEED experiments with various nitrogen coverages 
and beam energies. They found that a fractional coverage of approxi- 
mately 0.4 (defined as the ratio of the adatoms to the number of sub- 
Strate atoms for a specific face of the target) was relatively stable. 
At higher coverages, of about 0.5, the nitrogen would be absorbed into 
the bulk at temperatures of about 1000 XK. The 0.4 coverages were 
Stable at this temperature. The model suggested by Griffiths, et. al. 
was a contracted domain structure, where the nitrogen occupies a four- 
fold hollow site, but that the four surface tungsten atoms are uni- 
formly displaced towards the adatom. This leads to a seres of "is- 
lands" consisting of 16 nitrogen atoms and the corresponding surface 
tungsten atoms, as shown in Figure 4. This type of formation was phys- 
ically modeled by conducting a series of laser diffraction experiments. 


Laser diffraction gratings were constructed, with an average island 


Zo 


size Of 4x4, and when tlluminated by taser light produced diffraction 
patterns similar to those of the LEED beams on the actual target. 
Again, while this provided insight tnto the relative placement of the 
nitrogen on the surface, it did not provide any information on how the 


nitrogen was related vertically to the crystal surface. 





llustrates contracted-domain structure of the 0.4 
monolayer of Nitrogen on W(i00) or Mo(100) surface. 
Large hatched circles illustrate the top layer of the 
substrate atoms, the small filled circles, the 
mtrogen atoms. 

a) Plan view: Shows Domain and Boundary of Structure. 
b) Cross Section through line AA. {[Ref. 54, p.lis0em 


Figure 4, Suggested "island" Formation of Nitrogen on Tungsten. 
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Meyernoff assumed the equilbrium position of the nitrogen in his 
researcn. Webb examined equilibrium points both above and below the 
Substrate surface plane. While there was a great body of information 
available on the adsorption of nitrogen on single crystal surfaces, 
there was little data on the sputtering of nitrogen from. single crystal 


Surtaces. This was rectified by a recent report by Winters [Ref. 55]. 


B. WINTERS’ SINGLE CRYSTAL SPUTTERING EXPERIMENTS 

This paper (Ref. 55] was a2 continuation of Winters’ previous 
Studies in the sputtering of chemisorbed nitrogen on tungsten. In fs 
1974 paper he proposed two mechanisms for the sputtering of nitrogen 
from the tungsten surface (Ref. 27]. He suggested that for low energy 
incident wns, the primary mechanism for the sputtering of the nitrogen 
was direct Knock-on collisions with the impinging or reflected bombard- 
ING IONS. The proposed mechanism for sputtering at higher bombarding 
energies included nitrogen atoms Knocked away by sputtered substrate 
atoms. Winters’ 1982 paper ([Ref. 48} examined the sputtering of 
nitrogen from molybdenum and tungsten polycrystalline targets. 

This latest report by Winters and Taglauer described experiments 
investigating the sputtering of chemisorbed nitrogen from  W(i00), 
W(111), W410) and Mo(i00) single crystal surfaces. The targets were 
bombarded with helium, argon and xenon ions in the energy range of 300 
ev to 5000 eV. The relation between the sorbate mass, substrate mass 
and sputter yield were examined. Conclusions were drawn about thermal 
spikes, recoil implantation, and cascade = mixing. since these exper- 


iments were conducted on single crystal surfaces they provided some 
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basic mnformation on the physics of the sputtering of muiticomponent 
systems. Additionally, since the experimental data was derived from 
single crystals, they lend themselves more directly to computer model- 
Ing. A bnmef description of the experiment indicates how the system 
can be modeled in the computer simulation. 

A Faraday cup and the four single crystal targets, W(100), W(‘11), 
w(110), and Mo(i00) were mounted on a rotating manipulator. The tar- 
gets were cleaned by heating to about 2500K for tungsten and 2200K for 
molyodenum by RF induction heating [Ref. 56}. The sample was exposed 
to ION» until about 1/2 monolayer was adsorbed. This iS in agreement 
with the findings of stable nitrogen adsorption on tungsten as men- 
tioned in references 5e thru 54. The use of the ‘Nn isotope allowed 
one to isolate the adsorbed species involved in the study. 

After the adsorption, an Auger spectrum was taken to ensure the 
absence of impurities from the sample surface. The ion beam was 
scanned by the Faraday cup in order to ensure beam uniformity, and to 
measure the beam current density. The ion beam was rastered across a 
0.48 cm aperture, impinging at normal incidence upon the sample cen- 
tered behind the aperture. A mitrogen Auger peak-to-peak intensity was 
determined as a function of the ion dose. Finally the sample, which 
was partially cleaned, was exposed to the ambient for approximately 600 
sec after each run to monitor the readsorption of nitrogen. If read- 
sorption was observed, the run was discarded. 

The care taken in these experiments is quite apparent. The 
parameters of this experiment lend themselves quite well to computer 


simulation. Winters’ also presented a comparison of the experimental 
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results with theoretical calculations based on the Sigmund-Winters 
theory for the sputtering of chemisorbed gas [Ref. 27}. One recommend- 
ation in the paper was 

In particular, it would be useful to compare an extensive set of 


molecular dynamics calculations with both the experimental data and 
also the calculations presented in this paper. [Ref. 55, p. 26} 


C. THESIS OBJECTIVES 

The primary obyective of this thesis 18 to examine sputtering from 
the (100) surfaces of single crystals of molybdenum and tungsten. 
Nitrogen will be placed in different equilibrium positions in-~ relation 
to the substrate surface, and the system will be bombarded with argon 
ions of various energies at normal incidence. The effect of the mass 


Of the substrate in relation to the sputtering yield will be examined. 


La 


1t1. COMPUTER SIMULATIONSAND MODE DL DEN Siete 


We THE COMPUTER MODEL AND RELATED PROGRAMS 
1. QDYN86 

The computer simulation used in this study is called QDYN66, 
which stands for the 1986 revision of the QDYN (Quick DYNamics) program 
used by Harrison at the Naval Postgraduate School. This program uses 
multiple-interaction (MI) ftogic in & time-step approach. ine TU imeeae 
logic 1S appropriate when many events occur simultaneously. The model 
iS based on classical mechanics, using Newton’s laws expressed in Hamil- 
tonian form to simplify the calculations. The initial inputs to the 
program include the target crystalline structure, atomic messes and 
potential functions of the substrate, incident ion, and any adatoms; 
adatom locations, and bombarding ion angle of incidence, energy and 
impact point. The program develops the subsequent collision cascade, 
tracking the positions and velocities of the target atoms through time. 
The computations terminate when there 1s insufficient energy for any 
further ejections to occur. 

The positions and velocities of the moving atoms are tracked 
through each time-step. In order to strike a compromise between exces- 
sive computer run time, and yet to maintain reasonable energy conserva- 
tion, the time-step increment is variable. The time-step 15 determined 


by a specified distance divided by the highest atomic velocity in the 
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target. ive Orso tonce Chosen 15 0.1 Iattice units (LU), where the lat- 
Prceasintt i> GCeTIned ais one-Nalr the lattice parameter, a, (for a cubie 
Siructure). The time-step increment is further modified by other fac- 
tors which take into consideration the previous velocity, smoothing the 
transition in time. 

The velocities and positions of every atom in the target are 
not calculated for each atom at every time-step. Again, this would 
consume excessive amounts of computer time. As a result, the forces are 
not computed on a specific target atom until the atom «+s struck by a 
moving atom. Atoms which rise above the target surface are put in a 
tentative list of ejected atoms. After further atom ejgection from the 
target is unlikely, the atoms in this J18t are tested to see if they 
were travelling with sufficient velocity to overcome the attractive 
forces of the target atoms. Atoms which fail this test are added to the 
target, and are not counted as having been sputtered. Specific actions 
are taken for atoms that leave the bottom of the target, and through the 
sides. The program maintains, and per:odically updates a listing of an 
atoms nearest neighbors. After the completion of a collision cascade, 
the program reinitializes, and prepares for another trajectory. A num- 
Der of trayectories are run at different impact points in order to 
obtain better statistical results. 

Subsequent sections in thts chapter will provide further 
details on the choice of potential functions, substrate characteristics, 
impact pornts, and other factors as they pertain to this study. A num- 
Der of other articles provide more detailed information on the QDYN 


Simulation, and other molecular dynamics simulations. Reference 36 


Ze 


provides a good overview of sputtering models. Harrison and Jakas com- 
pared the MI and binary collision models fRef. 57}, and with Webb, 
examined a hybrid code between the two models [{Ref. 58}. The most 
recent, and most detailed discussion of the QDYN program by Harrison and 
Jakas identified uses of the program for studies involving insulators 
and semiconductors [{Ref. 59]. 
Ce. Ancillary Programs 
A number of other programs are used to support the main simu- 
lation. These programs are identified and are briefly described below. 
a. QDYNLIBB (QD86L! BB) 

This program forms a jibrary from which the main program 
derives routines for calculating the forces and potentials. The poten- 
tial and forces tables are dimensioned to 10000, and support four poten- 
tial functions. Adatom routines are available for placing adatoms on 
the surface of the target, or creating steps or vacancies. 

b. TARGLIBB 

This program forms a library of routines used for generat- 
ing (100), (110), and (111) faces of BCC, FCC and Diamond lattice tar- 
gets. A rotated FCC (001) lattice can also be generated. 

c. ANMOL 

This program analyzes the basic data set generated by the 
main program. The output of the main program consists mainly of a lJist- 
ing of the positions and velocities of the atoms. ANMOL uses this data 
file to calculate which atoms are sputtered, from which layer in the 


target lattice they originated, and rotational and vibrational analysis 
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of the ejyected atoms. Determination of multimers, ejection time dis 
tributions and eyection energres are also calculated. 
d. ANPLOT 
This «+S 38 program used toa interface with the DISSPLA 
Qraphics package. ANPLOT is used to generate graphical analyses of the 
data file generated by the main program. 
e. POTTER 
Thtis program is used to generate potential functions. An 
abbreviated input deck is read by POTTER, and two data files are gener- 
ated. These fries, FORCE DATA and ENERGY DATA, are used by the next 
program to obtain a graphical representation of the functions. 
f Se Owr Om 
This «tS another graphics routine that interfaces with the 
DISSPLA graphics package. POTPLOT uses the data files generated by 
POTTER, and produces a liner plot of the potentials, and a semi-!log plot 
of the forces. This allows one to quickly determine which modifications 
must be made to obtain the proper shaped curve for the desired potential 
and force function. The potential functions used in this thesis were 


generated using this program. 


ae SUBSTRATE AND ADATOM PROPERTIES 

The systems of nitrogen-on-molybdenum and nitrogen-on-tungsten were 
chosen by Winters for his sputtering experiments for avery specirfic 
reason. Many of the physical parameters of the two substrate materials 
are very similar, and the behavior of nitrogen with each was thought to 


be almost identical also. This provided a means of physically reducing 
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the variables in the system, enabling him to s:so0late and concentrate on 
the mass difference of the two materials. A summery of the Key physical 


parameters for the two systems 1s show below in Table 1. 


TABLE 41. PHYSICAL DATA FOR MOLYBDENUM AND TUNGSTEN. 


MOLYBDENUM TUNGSTEN 
Atomic Number (Z) 42 14 
Atomic Weight (amu) 95.94 183265 
Atomic Radius (A) i503 187 4 
Density (gm/crr) 10.22 19.3 
Crystal! Type BEG BCC 
Lattice Constant (A) 3. aT 3.165 
Lattice Unie=(17eadageen tools 1.95825 
Re (A) e.8 2.894 
Valence +4 +4 
fonic Radius (A) O70 Oe 7Ae 
Cohesive Energy (eV) 6.82 8.90 
Binding Energy of N (eV) B35) 6:5 


Data for Table 1 derived from References 60-62. 


Table 1 illustrates the high degree of similarity between some of 
the physical properties of molybdenum and tungsten. Nitrogen is assumed 
to occupy the four-fold position an the surface of both materials. The 
major difference between the two materials is the significant difference 
in’ mass and =<Z, The pertinent physical properties of the adatom, 


nitrogen; and the incident ion, argon; are summarized below in Table e. 


TABLE 2. PHYSICAL DATA FOR NITROGEN AND ARGON. 


Nitrogen Argon 
Atomic Number (Z) i 18 
Atomic Weight (amu) eo) 39.94 


Data for Table e@ derived from Reference 63. 
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om POTENTIAL FUNCTIONS 
‘ia Gener 4a | 
The dynamics of the atomic interactions are controlled by the 
enoice of the potential functions used to describe’ the potential 
Breer Oies and forces felt Dy the atoms. There is no "right choice" for a 
potential function in these types of simulations, indeed, Harrison nas 
called the development of useful potential functions for simulations a 
"Dlack art". A number of different potential functions can and have 
been used in sputtering simulations. Harrison has found that the choice 
of potential function parameters do not overly affect the outcome of the 
Simulation [Ref. 64}. The potential functions used in this investiga- 
tion are described briefly below. Detailed discussions of these, and 
Gmmer potential functions can be found in Torrens [Ref. 65}. A number 
Of survey articles by Harrison describe more fully how the potential 
functions are applied to the simulations [Refs. 36, 57-59]. 
a. Born-Mayer 
The Born-Mayer potential function is a purely repulsive 
function. This function is used for intermediate atomic separations. 


The form af the function 15s 


V(r) = a exp(-br). ame 
b. Moliere 
The Moliere potential function 18 another purely repulsive 
DoeemtidhweratrcGt10n, Of a slightly different form. It 18 called a 
"Sereened Coulomb Potential," and is an approximation of the Thomas- 


Fermi screening function. The general form of the potential is 


DE 


V(r) = (24525 ec / r) {0.35 exp (-O.3r/a) + U. SogexXDieieen oye 


O.4 exp (-6.0 r/fa) J . ( cm 


The term "a" 15 called the "Firsov length" defined as follows: 


a = K 0.8853 ap /(Z, 174 + Zo '/ej)e/3 | (989 


The “ap" term is the Bohr radius, 


ap = he / 41& mee = 0.5292 A. ( as 


The "at parameter may be modified, since the Moliere potential is an 


approximation to the Thomas-Fermi function. When k = 1, the potential 


iS pure Moliere; when K is set unequal to 1 then the function 1s Known 


as a "modified Moliere" potential function. 


GC. _ =norse 
The Morse potential is both attractive and repulsive, 
depending upon the separation distance, "r". It has the form: 
V(r) = De exp [-e a (C= 6) == 2a exe C SaOG( Pisatis ) Jue ( om 
The term, De is the well depth, and fe 1s the equilibrium separation of 


an atomic pair, and alpha is a scale factor which controls the shape of 


D4 


ero eve ec TUCK Oy) 1S attractive for r > fe,7and repul- 
serve Tor rT < Po: 
d. Composite Morse-Moliere 

The composite Morse-Moliere potential function is used to 
more closely model the dynamics of the sputtering events. The repulsive 
wall of the Morse potential is jyotned with 3 cubic spline to the Molrere 
potential, to form a single potential function which can be used over 
long ranges. The two functions are joined by varying the alpha term of 
the Morse, and the "a" term of the Moliere to obtain an intersection 
which has a smooth, continuous slope. The repulsive wall governs the 


colliston dynamics, and the dttractive well of the Moliere controls the 


sputtering by determining whether an atom will escape the surface of the 
target. 
ae Selection of Potential Function Pdarameters 


a. Subdstrate-Substrate, Adatom-Adatom Function Parameters 

A compos:rte Morse-Moliere functton was used for the solid 
phase Mo-Mo and W-W potential functions. The specific parameters used 
in this thesis follow those developed by Webb, with one exception. The 
function used for the W-W forces had a discontinuity as shown in Figure 
oe This was smoothed by changing the Rg value. The potential curves 
are shown tn Figure 6. 

The nitrogen-nitrogen potential function is a pure Morse 
potential function. The tabulated data for the solid phase W-W, Mo-Mo, 


and N-N potential functions is listed below in Table 3. 
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TABLE 3. SOLID PHASE POTENTIAL PARADE 


e 
(eV) (A) (A!) (LU) (LU) (LU) 
Mo-Mo 0.997 2.800 1.519 ©. 0 0.790 0.836 2.40 
W—W RSS iS. 2.894 1.200 O20 0.80 121360 C wae 
N-N Wee StS 1109'S 2. f00 0.0 0.0 0.0 1 ea 


The a and kK terms were used to metch the slopes of the 
Morse and Moliere potential functions, as mentioned earlier. The value 
of Dewwas selecteo if Of Gen ate op cain the proper cohesive energy for the 
substrate as shown in Table 1. The Re term is the nearest neighbor 
distance. The spline boundaries, R,z and Rp, are in turn a function of 
the a and K values selected, and should be no more than ¥% LU apart. This 
topic will be covered in more detail in the following section. The 
final value, Re, iS a distance over which the force and potential 
calculations are made, a cut-off distance. This allows for tnteractions 
only between nearest and next-nearest neighbors. This value ts set at 
1.71 LU for all interactions except the substrate-substrate interac- 
tions, where it ts set at 2.40 LU. 

b. MO-N, W-N Potential Function Parameters 

Six new potential functions were generated to examine the 
effect of adatom placement on the sputtering cross section. The FORTRAN 
routines POTTER and POTPLOT were used to help generate these functions. 
The first step was to select an elevation to place the adatom, and from 
that to determine Ro, the nearest nelgnpor separation distance. Re 15 
the distance measured vertically between the cubic center and the ada- 


tom, as shown in Figure 7. These values were then s:nput to POTTER. 
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POTTER generates two data files used by the POTLOT rout:ne 
ap OC ie eePOUSiinalsana TOheces curves. fhe R,, Rp, K, and alpha 
mmmes ore modtiTIeCGd SO that 4 $mootTn curve is plotted for hoth the 
potentials and forces. Once this rough determination is completed, then 
Peelo meditied and evaluated by running a single trajectory at one 
time-step under CMS in order to obtain the proper sublimation energies. 
Pecemodis1cd sO thot the sublimation energy of the nitrogen to the 
Substrate 15 maintained at 6.5 eV. The binding energy of the adatom is 
one-half the sublimation energy, or 3.25 eV. 

After the sublimation energy has been properly set, then 
the revised values for De are once again run through POTTER and POTPLOT 
in order to insure that the curves have the necessary smooth = slopes. 
These graphical methods for rapidly evaluating potential function par- 
ameters have greatly eased the process for generating potential § func- 
tions. The values for the Mo-N and W-N potential functions used in this 


thesis gre ltisted below in Tables 4 and 5. 


Adatom 
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= Placement Height 
Substrate 


Figure 7. Physical Dimensions Used in Potential Function 
Calculations. 
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TABLE 4. Mo-N POTENTIAL FUNCTION PARAMETERS 











Location De Re a , K Rs Rb Ro 
(A) (eV) (A) (A) ~ (LU) (LU) (LU) 
~Q.05 2.970 1.524 2.60 0.80 0.540 0.560 1. 7 ieee 
C145 e2e4e5 1.720 eae) 0.80 0.540 0.560 {7 
Orcs. 1,997 12736 Cais 0.80 0.500 G2oce 1.71 
GreI ee te 1.864 2.20 @,.0 0.560 DO -oc8 1; 
G4555 als (ee 1.908 emcee O20 Oe oe 0.560 ‘. aa 
O23 1.60 Ee ele (6. 2.15 0.0 O7 reo 0.750 1.71 ** 
Note: Rows marked with ** indicate parameters developed Dy Webb, 
and used in both theses. Unmarked rows indicate new potential func- 
tion parameters. 
TABLE 5. W-N POTENTIAL FUNCTION PARAMETERS 
Location De Re aX K Ra Rb Re 
(A) (eV) (A) ee CEUs ae (LU) 
205 3.292 1.524 AaB) 0.8 0.65 O67 1.71 
0. tes 2.341 1. (446 e.te Oe, ORS G7E8 1.7m 
0.2464 22O\s 1.8289 2.60 Ore Q.80 O26 1. % 
©. ses 1a 64 1.9075 Cua Ore Orso 0.90 1.71 
O2467 ie 6s 2 1.960 e.. Se 0.9 O83 0.85 1.71 ex 
Note: Rows marked with x** indicate parameters developed by Webb, 
and used in both theses. Unmarked rows indicate new potential func- 


tion parameters. 
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C, lon-Substrate, !on-Adatom Functions. The potential func- 
tion parameters for the incident ion must be considered. The target is 
bombarded by Argon ions, Dut for the purposes of these calculations, the 
charge on the sion is neglected. The ton is charged so that it can be 
accelerated to the desired energy to bombard the target. Since the 
incident ton is a noble gas, the event that it will react with the 
molybdenum, tungsten or nitrogen atoms other that through simple colli- 
$10n process is unlikely. The potential functions are correspondingly 


modified Moliere. The potential function parameters are listed below in 


Table 6, 
TABLE 6. AR-N, AR-MO, AND AR-W POTENTIAL PARAMETERS 
(eV) (A) (A) 7! (LU) (LU) (LU) 
Ar-N 0.0 Ome ORO 0.0 ae Wee ree S| 
Ar-Mo 6) ©@ oO .0 OF 0.0 Veet cern seen e | 
Ar—-W 0.0 0.0 Ono 0.0 deer | 1. 74 decd 


Figure 8 shows the Nitrogen-Nitrogen potential function, and 
Figure 9 shows the Ar-N, Ar-Mo, and Ar-W potential functions. 
d. Vacuum Phase Potentials 
Once atoms have sputtered from the surface of the target, 
there is a probability that further interactions can occur in the vacuum 
phase. Vacuum phase interactions were not emphasized in this study, but 
for the sake of completeness, the vacuum phase parameters derived in 


Webb’s thesis are included itn Table 7. There is little available data 
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Figure 9. Ar-Mo, Ar-N, and AR-W interatomic Potential Functions. 
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for the Mo-Mo, W-W, MO-N, and W-N systems. The following substitutions 
were made, using similar systems to those examined in this study. The 
Ag-Br system was used fro Mo-Mo, and VO was used for Mo-N. Data on the 


N-N system was available and used. [Refs. 50, 65). 


TABLE 7. VACUUM PHASE POTENTIAL FUNCTION PARAMETERS 


DA a Re Q Qx 

(eV) (Ayes (A) 
Mo-Mo {. f00 1.790 2.390 247.000 0.6079 
N-N 5. 760 3.359 1.098 C305 .570 14.320 
Mo-N 6.100 1.965 1-324 1020.000 4.700 
W—-W 3,100 1.450 2.470 308.399 0.960 
W-N 4.900 2.159 ee Tae ia' 967.000 4650 
RHC 2.900 1.560 2.390 1050. 000 5.000 *x 


Note: Tne Rhodium vacuum data was used in some cases. 


D. TARGET and |MPACT AREA CONCERNS 
One characteristic of MI Simulations of this type, 1s that finite 
target s:zes are used. Two concerns must be balanced in determining the 


optimum target size. First one must try to have a target of sufficient 


size such that the majority of the sputtering events are contained. lhe 
the target :s too small, then information will be lost due to a failure 
Of COntamment Seven. On the other hand, an excessively large target 
will require excessive computer time to complete each trayectory. In 


this thesis, two target sizes were used. At the SOO eV energy level, 2 
19x6x19 target was used, and at higher energy levels a 23x8xe3 target 
size was used. Webb used the larger target in his thesis. 

Nitrogen atoms were placed in the four-fold position on both target 


sizes. 40 atoms were placed on the smaller target, and 60 atoms were 


Ad 


Biacea ‘on the larger. In both cases this resulted ina fractional cov- 
erage of approximately 0.4, which corresponds to the contracted domain 
eiructure of Reference 94. The nitrogen adatom coverage is shown in 
Figure 10, corresponding to the (001) face of the substrate material. 
Table 8 tists target area and nitrogen coverage parameters. 


TABLE 8. TARGET AREA AND ADATOM COVERAGE PARAMETERS 


Target Size 19 x 6 x 19 ao xX oex 23 
Substrate Mo W Mo W 
No. Nitrogen 40 40 60 60 
Area (cm) 8.02x107'4 8.11x107'4 1,.20x107'13 4,24x10713 
Fractional 

Coverage oe. 0.4 0.42 0.42 
Coverage 

(atoms /cm®) 4.99x1014 4.93x1014 5.0x10'4 4.96x1014 


Another factor which impacts on the run time of the simulation, and 
the statistical accuracy of the results is the number of impact points, 
or trajectories run. The impact points are designated by a pre-selected 
"Din-point", and then the trajectories impact in a representative area 
based on this point. Two types of pin-points were used. in one, the 
"nit" case, a nitrogen adatom is placed in the representative area, and 
can be hit directly by the incoming ion. The other, "“non-hit" case, 
does not have a nitrogen atom in the representative area. The simula- 
tions were run in this study with 150 trajectories per pin-point, or 300 
trajyectories per complete run. Previous studies used 300 trajectories 
per pin-point, for a total of 600 trayectories per case. 

The pin-points for the "hit" case were located at (8.000, 6.000) 


for the 19x6x19 target lattice, and (13.000, 15.000) for the 23x8x2é2 
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target. The non-hit cases were run with pin-points at (8.000, 8.000) 


and (11.000, 11.000) respectively. These pin-points are measured on the 


K= 


and z- axes. The representative areas and pin-points are also Shown 
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smal] circles represent nitrogen atoms. 
"H" indicates "Hit" representative area. 
"NH" jndicates non-hit representative area. 


Figure 10. Top View of Nitrogen Reacted Surrace. 
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eee micta iim) tlie | ON 
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Me eee Wee 26C T PONS 
Pye] comeept ar the sputtering cross See antl Te ey To the 
aria lyses im this study, ine deryvyatiOnentT the cross Section is carried 


Out in detail in Appendix A, The basic methodology used by Winters was 


adapted to the simulation, and the following relation was derived: 


Thy - - fin (1 | Y/No) 3 A, ' { a ) 


eee vue tet  yrclad thom the simlation, oy to the Nitrogen sputter - 
ae Creer oe eo C1 OFM, No 6 the initial number of adatoms placed on the 
target and Ae Leawvicerunrroce area tie target, 
en witeers Res u 1s 

Winters? experiment involved the bombardment af the (100) 
Plane of Molybdenum and tungsten with argon, xenon and helium pons mn 
energies from 300 eV te 5 KeV. The caszes Studied in this thesis con- 
Pete argon pons Of the SOU €V To ¢€ KeV Fangde of MO(100) and w(10Q). 
Wanecer S? cre ctimere fesud lis 11 om Nis SxpDemmnents, and tireoretical calcu- 
lations from the Sigmund-Winters theory of multicomponent sputtering. 
fewer toe that fall wri the botnds OF ths stucyucare tabulated 


fevewerr lace 2 (RET, 5S p. e9). 


4/ 


TABLE 9. WINTERS’ THEORETICAL BND EXPER IMENTAL Osa ce oer 
AS A FUNCTION @F [ON Sar hia) 


hon Exper iment. Calculation 
Energy (1074® cme} (1071 eme) 

(eV) Mo (100) W100) Mo(tQO0) W(10Q0) 
900 are e AG &.O r.9 
1000 6.5 8.6 &.8 T.4 
COO0O A,B 9.0 8.0 Reset 
3000 AG &,6 &.9 B.6 

3, Previous Theses at the Naval Postaraduate School 


Sy Meverhotrt 

Meyerhoff's thes1s «1 1983 examined clean and hitrogen 
reacted molybdenum and tungsten (001) surfaces bombarded with argon s0ns 
of energies ranging from 500 &V to 3 keV. WNetrogéen was placed “ie 
four-fold position, and was @xamined at two elevations above the (001) 
surfaces, The effect of the mass of the substrate was examined by 
replacing the Mo-mass with the tungsten mass. Meyerhoff concluded that 
the sputtering cross section of the nitrogen was more dependant upon the 
distance of the adatoms above the surface then upon the mass of The sub- 


strate. (Ref. 49] 


Webb, if early 1986, continued the study of fitrogen on 
molybdenum and tungsten. He refined the potential functions used by 
Meyerhoff, and performed basic sputtering research using the QDYN Simu- 
lat rons He studied the sputtering cross section of nitrogen for place- 


ment at three different heights in relation to The (001) Surfaces im 
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Molytdenum yad Tungsten. eee ip ta! Lhe rip lregqen was me tf fibe-ly 
Loca LS e ewe toate, S25 °0NPdsed To heing ih the well formed by the 
mop tace atoms. 


4. Analyses Conducted in This Stud 





Tie Seep emeat ere oghalyses conmogucted in This study differed 
somewhat from previous theses. The previous studies did not differen- 
recom het wee te etree ts OF THE WHeiGent ran intially Hitting the 
geeorbed nitroden, “er the substrate. The Simulations were run with 
miter pots Coverrnd beth cases, but The results were not looked aT 
Berar ately. ities seedy, Bhe tTWe cases wwere analyzed separately. 
tee was «One in Order fo gain Insight into possible sputtering mech3- 
eos aga OUTHINEed earlier. 

Daw ocener Leland Cross esee tions Were examined at 300 eV 
Greener ay sputtermdad) afd ec keV (high enerdy sputtering). ive 1) 1 1 eat 
So teen Weleeeaanined in Order [Lo determine the sputtering mecha- 
nism of the mitrogen. The yield of the substrate was examined, and the 


Paomoretiageparme Lo the reacted suDSTmMate yrelWds Was Used 45 4 measure 


CIs 


meeerenrning tie relative effects of the sdsorbed nitrogen. The purpose 
eens SGI YS1S Was to déetermine if the mitreadgen enhanced the sputter- 
foe Geer the NMrtrogen, particularly at the higher ion energies. 

he Cthect Ofeatne Giss Of The substrate was examined, in order 
tO gttempt to resolve the difference between the findings of Winters and 
Webo and Meyerhoff. The technique used by Méeyerhoff was again used in 
oe ht tidak the sOULTEr ind croas sectron of ni troqenr™was 
evaluated at three adatom placements as opposed To one, and the hit, 


= = -— = 


no-hit cases were examined separately. 
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EB, MOL FEDENUM RESULTS 


Nitrogen Souttering Cross Sections and Adatom Placement. 





4, 200 é€V Jon Bombardment 

The sputtering of nitrogen 563 & function of adatom place- 
ment was examined at two energy levels, 500 &V and 2 keV, lr addi t ae 
the “hit" and “no-hit" cases were examined separately, Table 10 Jista 
The nitragen sputter yield, and cross section 48 38 function of adatom 
Placement for The 500 €V bombardment. The results are plotted in Figure 
114 for the fit, no-hit and Total cross sections. Tne plot of cross $6@e¢- 
Tion vs atom Placement results in aM apparent “well" corresponding to an 
adatom placement at 0.245 A above the (100) surface, This Was trueiiaa 
both the hit and total cross section cases. The sputtering Ye) Gee 
micreased 45 4a function of 4datom distance from the surface. The simula- 
Lion cross section at this point 18 5,63 x 10795 cm®, corresponding to 
the experimental value obtained by Winters of 5,6 % 10 7'9 cme. 

The effect of the ion Hitting The nitrogen 1s guite 
agfparent. The sputtering cross section of the hit case increase by 4 
factor of frome to 4 sbove The non-hit case, and by a factor of about 
{1.5 for the total cross section. This indicates that the primary Meche 
nism for the sputtering of the nitrogen 1% when the incident joie 
The adatom, and it in turn 13 reflected from the substrate, oS ies 
treated eat tie iy Figure 2(a). This seems to be more pronounced ™= ia 
adatom i3 placed further above the "well" at PCa Ae The adreemert 
between the experimental and simulation Ppesurts for’ this point “nde 
That the nitrogen may be located 3f approximately that height above The 


(O01) surface. The nitrogen yields and cross sections are listed Gbelew 
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Paaure. 11. 


ADATOM PLACEMENT (ANGSTROMS) 


Comparison of “Hit", “No-hit" and Total Cross Sections 
With Theory and Experimental Values for 50Q €V Ar on 
Molybdenum. 


lk 


with H corresponding to the “Nit “Case iNet Or Uilecueeeetel eee and f 
for the total values. 


TABLE 10. NITROGEN SPUTTE RMD TA, MO(O01) ove 


Adatom Spu Uter Cross 
Placement Yield Section 
(A) (x 107-16 cme) 
H 0.41 9.4 
-QO.05 N Oe he. 4.3 
T 0.28 6.7 
H OFoc. vO" D 
0.146 N Ocalks ec .68 
ol Oras OCT 
H QO.42 8.47 
OQ ..2a N 0.14 C..0u 
T 0.28 Soo 
H 0.63 — 12.8 
QO.290 N O20 4.02 
T Oe4e 8.4 
H 0.66 toiens 
O7s30 N Ove 4.02 
if 8.67 
H Ono 16.3 
O.380 N O.19 375 
ai QO.50 Hone 


Db. 2& KeV Bombardment 
Nitrogen sputtering when bombarded at 2 keV was examined 
at the same adatom locations, and the data 1s tabulated itn Table 11. 
Figure 12 shows the nitrogen sputtering cross sections as a function of 
adatom location, with the hit and no-hit case considered at four loca- 
tions. The curves follow the same general trends as Figure 11, with some 


exceptions. The same "well" is observed, but in this case, the _ total 


a2 
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LEGEND 
Oo TOTAL CROSS SEC 
Cet” CeCsae 
oO SsigeCmoss SEC 
CALC: B.0 E-16 





we 
A... e 
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ADATOM PLACEMENT (ANGSTROMS) 
COMpar ison of "Hit", "No-hit" and Tote} Cross Sections 


With Theory and Experimental Values for e K*V Ar on 
Molyodenum, 


De 


0.5 


cross section of 5.6 x 107'% cme 1s below the experimental and theor eam 
cal valWies=@e Sroe~ 10°18 cme and 8.0 x 10718 cre& respectively. The 
slope of the sputtering cross section curve seems to level out at adatom 
placements above 0.29 A. The cross section for the -O.05 placement 15 
below the “well” at 0.245 A. The -Q.05 point was plotted using Webb's 
data. Both he and Meyerhoff used 14 in their computations, where 13N 
was used in this study, and in the experiments by Winters. The 7Z mass 
difference may account for part of the difference in the result. The 
effect of the hit case iS again quite apparent, with significantly 
higher sputter cross sections resulting from the 10n StrikK 1) 


nitrogen atom. 


TABLE 11. MOEN) SPUTTERING CROSS) Chen TONGme Gmiuay 


Adatom Sputter Sputter 
Placement Yield Cross Section 
(A) (xX 10718 cme) 
-0.05 T 0.26 Ome 
H 0.56 (ieee 
0.146 N OTS S00) 
T O35 OH 
H 0.45 9.02 
OO; 245 N Ora, 2.2 
T Q.28 5:6 
H 0.79 15.9 
0.290 N 0.20 4.0 
T 0.50 10.03 
H 0.80 Out 
O7335 N 0.23 4.60 
i OFS ines 
0.380 il ORo3 10.6 
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Ca Poe Ween ete ST Ler ing Cross sect mans 


ie Svea sDULTLer ind GFrGSs sections for SOO «V and ¢& keV 


es | 


are plotted on the same graph in Figure 4 The general trends are 
apparent. Both curves exhibit the same "well" corresponding to an adatom 
placement of 0.245 A. The sputtering cross sections then increase as 264 
Function of iIncredsing herght of adatom placement, There seems to be 
something speéci4|l @bout the 0,2e45 A placement. Since the sitmlation 
Pee ose StGnOn agreed so closely with the experiment, +t 1S 4 good indi- 
meeronm that this mant be close to the actual location of the atom. The 


fact that the ¢e@ KeV curve 380 closely Matches the SOQ eV curve gives fur- 





ther credance to this assumotion, 
ae Nitrogen sputtering Cross sections a5 4 Function of Enerc 





The sputtering of nitrogen 48 8 function of energy was 
meamined. the purpose of this analys:s was To compare the results of 
Stes simlation to the data presented by Winters 38 listed in Table 9. 
Pee ecco Dlaitcdra: far energy levels, SOO eV and 1.0, 2.0, and 3.0 
key, The analysis was conducted for adatom placement at 0,245 A. This 
Pore WasecNOse, S1nGce if 560 Closely Matched the experimental! cross 
Section aS Mentioned above, The results are shown in Figure 14, with 
fee COrrespaondind experimental and theoretical values from Winters’ 
Paper, The values at 300 éV and 3.0 keV aré within the bounds of the 
Poteet ana theoretical values, but the points at 1 KeV and e keV 
Beeeeseyonaeeene limitsae The Variation of the simulation resalts from 
Pus OT Witten nniameiiedahe uno; the “Uncertainty of the siftwilatron 
model, However, the general trend of the simulation cross sections de 


Oia dal awit hewiiihers! ese s.. 
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NITROGEN SPUTTER CROSS SECTIONS 
5 AND 2 KEV AR<001> ON MO(001) 
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Figuré 13. Comparison of Total Cross sections for 500 «V and 2 keV Ar 
on Molybdenum, 
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jee eee NIALL RESULTS 


i j j 
0.0 0.5 1.0 1.5 2.0 2.5 3:0 aS 
ION ENERGY (KEV) 
Nitrogen Sputtering Cross Section a8 &@ Function of Energy 


for Adatoms Placed at 0,245 A on MO(I00), 
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Of The sous aie 

The ratio of the reacted substrate yield to the bare substrate 
yield was examined 47 SOO ev and 2@ keV. The hit, peat and tots) 
Vieldsmyere Consider eq, The results are listed in Vieabie Le, and are 
DTGT Lea ies eal: Ih alt cases, for a3) adatom placements examined, 
the yreld of the reacted surface was lower than that of the bare sur- 
face. Additionally, the substrate yrelds of the hit case were lower 
tan for the total snd non-hit case. This indicates thot the nitrogen 
adatoms absorb some of the energy of the incident particle, and there- 
fore Jess momentum 13 available to anitiate the collision cascade in the 
SUD Sia ane. 

The case Was the Samerat coKeyv- The sputtering yields OF eine 
Substrate were lower in é€ach case when the adatom was init~ally Sy ae 
Dy the ion. The incident 1on at 2 kev still failed to Nave Sit, ioe 
energy to impart enough momentum to the adatom for if to act 4% 8 Secord 
projgectileé and thereby enhance the sputter yield @F the S05 Gi (Cee 
yields of The substrate and the ratios with the baré substrate for 2@ kev 


are also listed in Table ite. 


= TUNGSTEN RESULTS 


pe Nitrogen Sputtering Cross Sections and Adatom Placement 





4. 500 é@V lon Bombardment 
The sputtering of nitrogen from (001) tungsten was 
examined at 500 eV and 2 keV, The hit and no-hit scenarios were 


examined as with molybdenum. The results ere listed in [a0)e tee 
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RATIO OF REACTED TO BARE YIELDS 
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Figure 15. Ratio of the sputter Yield of the Reactegd oubstrate to the 
Yield of the Gare Substrate. 
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ae si! eee The Same general pattern wae Poteeris. tier 
molybdenum. Mie ees ross section far the hit scenarie 15 again 
Sragnificantly greater than for the no-hit case, for all adatom place- 
ments. Therefore, a5 with molybdenum, the majority of the nitrogen is 
Peeeeoeeg ao 4 result of & Girfect Collisian by the imeident © ian. hie 
Lo ks | SpuTTEP ING cross sections generally follow the experimental val- 
Wes until the point corresponding to the 0.2464 A adatom location, where 
roms hope of the cross Section Plot rapidly increases. 
b. 2 keV lon Bombardment 

The sputtering cross section and yield data for the 2 keV 
bombardment 15 sunmmeri2z2éed in Table 14 and plotted in Figure 17. [ETS = 
rs no "well" aS was abserved for the molybdenum case, Nowever there is a 
feeanethn the slope of the curve at the point corresponding to adatom 
Placement at 0.2464 A. While not plotted, the sputtering cross sections 
Tar stne Nit case were aghin significantly greater in all instances than 


for the non-hit case. 


¢c, Comparison of Sputtering Cross Sections 
The total nitrogen sputtering cross sections for the 500 
éV and e@ KéeV bombardments are plotted in Figure 18. The two curves fol- 
few the safe general trends. aowe sit eitileenoeer ved, DULL there i165 an 
Siete enon Dpeint mm the curves at a Point corresponding to ah adatom 


Piacememeueor ©.¢464 A. One special Nete, the Spetteringd cross section 
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Nitrogen Sputtering Cross Section 45 6 Function of Adatom 
Placement for 2 KéV Aragon on W(100), 
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NITROGEN SPUTTER CROSS SECTIONS 
5 AND 2 KEV AR<O01> ON W(001) 
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Figure 18. Comparison of Total Nitrogen Sputtering Cross 
Sections for 500 eV and 2 keV Argon on W(100). 
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TABLE 13. Wit10O0) NITROGEN SPUTTERING FESULIS, SQQ =y 


Adatom abutter Sputter 

Placement Yreld Sore See 1OF 
(A) (x 10719 eme) 

H O.60 Teno 

— iD he. N 0,23 4.68 

T QO, 42 a als, 

H 0,69 14, 4 

O,1623 N O.e2 4,68 

T @e46 o) 5 Bie: 

H 0,68 i329 

00,2464 N Q0,27 5,49 

T O.47 9.909 

H 0,83 if ..© 

DAs ape) N 0,32 oe | 

Nem aot 11.6 

H 0,99 60.3 

Oa vows N Q.33 S272 

jae eG Noe a 


TABLE 14. NITROGEN YIELD AND SPUTTERING CROSS SECTIONS FROM TUNGSTEN 


2,0 keV 
Adatom Tota | SouiLiLer 
Placement sputter Crass section 
(A) Yield (x 10718 eme) 
-~.05 0,28 5.67 
.1623 0,35 7.08 
2464 0.36 7.26 
325 0.40 8.09 
487 0.45 ei 
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of mitragen aT 500 #8V ys S1gni Piet ly on eater thai ac ey Gl ee 


adatom placements. This +5 similar to the findings in the two pPrevigie 


theses [Ref. 49, S50]. This conflicts with the experimental anew thadrer 


ica} values presented by Winters. This 18 an indication that the mode} 


used in the simulation 1:8 not correctly portraying the behavior of the 


er. nee SPU TLE aS ane oie Sees 





in an effort to find an adatom placement Jlocation thal «agreed 


with the experimental results, the total nitrogen sputtering cro 


cn 
Cc 
co 
ia @ 
m 
| 


tion was plotted as a function of energy for 331 five adatom locations 
Lester The results are shown in Figure 19. in all cases, there Waa 
strong negative slope, that did not correspond well with the slight pos- 
itive slope of the experimental and theoretical values. Again thi Secor 
responds to Webb’s and Meyerhoff’s findings. We indicates Thal me 
potential! funeétLion My Net He (penteewiia the true behavior One 
Nitregery on The winds. cic 

a = of the 


mutter ing aubstrate 





The analysis of the ratvo of “they sputtering yietG Of aaa 
reacted substrate to the bare substrate was not conducted to the same 
détai) a5 with molybdenum. The results for he ¢asé3 examined are awe 
marized below in Table 45, In all cases, with the exception of the hit 
case for placement below the surface at 500 eV, the yield of the r é4 Gee 
Substrate was lower than that of the bare substrate. This was the case 
both at S00 eV and 2@ keV. This indicates that for these two levels, the 
adatoms again reduce the energy of the ion such that fess momentum 15 


transferred to the substrate, and the yield of the substrate i5 lowered. 
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NITROGEN SPUTTER CROSS SECTIONS 
AS A FUNCTION OF ENERGY, W(001) 
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Figure 19, Comparison of Nitrogen Sputtering Cross sections as 
a Function of Energy and Adatom Location, W(1i00). 
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One wnteresting observation 15 That The yi) ett Tie ep eee 
did inegrease, a8 expected, at The higher eneray, This despite The eae 
that the simulation cross section of nitrogen sputtering was higher at 
900 €V than at 2 keV. 


TABLE 15, SPUTTER YIELD OF SUBSTRATE W( 190) 


DON éV ec Kev 
Adatom 
Placement Yiélad Ynu/Yp Yoeld Yy/Ype 
A - 
ty 2,67 
Bare N Seo NA 

di c.98 2.63 NA 
id 1 ere 0,67 c.16 0.82 

-0,05 N etl 0,53 Seam f 0,60 
hi eae ie 0,59 1 aa QO,71 
H 4,83 0,68 

0, 123 N 12G6 Osa: 
ih 1,85 0.62 e.06 0,78 
HY eT We cu. 

0,244 N {,9 Gags 
T 1,84 0,62 Cue OF ee 
a Lee 0,66 

Seo N 2.05 O23 
T 1.90 0,64 eT oe 4a 
a ivi QO,65 
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we COMPARISON OF MOLYBDENUM AND TUNGSTEN RESULTS 


Comparison of Cross Sections 





The sputtering cross sections of nitrogen from molybdenum and 
tungsten are plotted in Figure 20, with the adatom placement normalized 


LO Fat Lice s . The molybdenum curves exhibit the "well" at 0,156 LU, 
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Figure 20. Comparison of Nitrogen Sputtering Cross Sections 
for Mo(100) and W(100) at 500 eV and 2 keV, Adatom 
Placement Normalized to Lattice Units. 


Ge 


but the fungsten curves Show a er iti 1 i ee yy ey 
relative point. ThrsS Point corpeepon@e 160 OF esa Ae) »bybdenum, and 
0.2464 A for tungsten. The potential function for the 0.2464. Een 
ment was dérived so That the behavior at The point corresponding to the 
molybdenum well could be examined. Again the major difference between 
the two substrates is the fact that the Sputterig cross ee1 oe 
ev 16 Higher than at 2 Kev for Tunes ten. 

The Hrt and no-hit dnalysts was conducted for Tungs (ene 
and the samé general trend was noted, Significantly higher sii eee 
cross Sections were noted for the hit scenario. ThtS Was Drie 
energy levels and adatom placements examined, THiS indicates tha 
primary mechanism for the sputtering of nitrogen if due to GiPect aie 
Suwon by The 1NGi Gem Ta, 

on Of SUDS Ua lent ne nme ale aes 


Compar 1so0n 





The ratios of the sputtering yrelids ofethesreseted Sua 
to the bare substrates at S00 €V and 2 keV bombardments indicate That 
the nitrogen does not enhance the sputtering of The substrate. This 15 
observed even at the relatively high energy level of 2 Bev. Both 
molybdenum and tungsten exhiorited this behavior. The incident i909, even 
at ec keV failed to transfer enough momentum to the nitrogen for +f to 


act 48 a second projectile to enhance the sputtering of The SUS) ae 


ee. DETERMINATION OF POSSIBLE MASS EFFECTS 
One of the major conclusions in Winters’? earlier works was the fact 
that the sputtering cross section seemed to be greater for the substrate 


with the higher mass [Ref, 48, 55), Meéyerhoff and Webb examined This 
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=) bee and concluded That adatorn placement had 4 greater effect. The 
fas a Teevewdamex canine if) This thesis by Selecting three adatom loca- 
mons Ii ene vreni yw tie Well aft O,c45> A, for the Mo(i00) surface, 
and replacing the mass of the molybdenum with the tungsten mass, leaving 
ee Our moO LS OF Ulimepotential fungtreans the same. The results are 
SlGiteq in Figure ei, for the fit, noe-hit and tote! cross section cases. 
eee Cgoc as, ce OULL er Wig Cr Gro seClhians were higher for the higher 
Peres MNGiCatco Uldteme hss OT tne Ssupstrate does have af effect 
Mme setELet id cross Section of the nitrogen. 

merece — siews the total nitragen sputtering cross sections for 
gee y MO(100), MOx(100) (The high-mass molybdenum), and W(100). TS 
araph 16 normalized for adatom placements aT approximately corresponding 
jocations. The trend for higher sputter cross sections FOr ew Sibel 
Poon oVious. Tivewieavy Moe has higner nitrogen sputtering cross sec- 
memes the true Molyogenum, and the tungsten hes the highest cross 
sections of the three. Through the use of the computer simulation, the 
feet feo ts Nave Heeheable to be truly rmolated, More so than wr pos- 
Sible ifn 4 physical system. This 15 one valuable aspect of Simulations. 

One poss1ble reason for the higher sputtering cross section for the 
feeier substrate ws that the dreater mass of the substrate provides 3 
better "springboard" for the nitrogen atoms. When the mitrogen 4aforns 
aieerneecrea trom their Pas+1tiens by the incident i49ns, they relatively 
hight nitrogen atoms will rebound with more initial momentum conserved 
from the heavier tungsten atams than from the Jighter molybdenum atoms. 

This model for the increase in sputtering yield with higher SLD 


Strate mass follows if the primary mechanism for the sputtering of the 
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Figure 21. Comparison of Nitrogen Sputtering Yields for 
MO(100) and Mox(100), 500 eV Argon. 
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Aetraden is direct CGblisiQn Gy” Vinee ee lf the mitrogen 4toms 
were Knocked off by reflected argon :0ns, or Dy SHhuttering substrate 
materials to 4&4 high degree, the effect of the mass of The substrate 


would seem to be less important. 


eS COMPARISON WITH PREVIOUS SIMULATIONS 

The results of this study indicates that therg 15 & 5) Ot) eee 
mass effect, this 15 contrary to thes fimdwngs of Webb “aid Shey eaiiaae 
The earlier studr1es did not isolate the effects of the direct ieee 
tion of the ion with the adatom to the degree that it was done in’ this 
& Tay The trend for the higher nitrogen sputtering cross section 4T 
the lower energy for tungsten 18 consistent through the three studies. 

The other aspects of the computer simulation, such 35 the angular 
and energy distribution plots remained consistent with earlier works. 
For g@ More detailed analysis of the “traditional” sputtering simulation 


analysis, the reader 1s referred to Repent eae iiaee2, 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The results of this simutation confirm the findings of Webb and 
Meyerhoff that the placement of the adatoms is most likely above’ the 
crystal surface. A location of approximately 0.245 A for the molybdenum 
Surface results in a nitrogen sputtering cross section that agrees with 
the experimental results of Winters. 

Comparison of the "hit" and "no-hit" scenarios indicates that the 
primary mechanism for the sputtering of the adatom is due to direct 
collisions with the incident ion. In all cases examined, the sputtering 
cross sections derived when the adatom was struck by the incident’ ion 
were larger than in the cases when they were not struck. This was found 
at both energy levels examined, 500 eV and e keV; and for both Mo(100) 
and W(100Q). This agrees partially with Winter’s conclusions, but does 
not show a significant increase in non-hit sputtering at higher energies 
he proposed. 

The mass effect of the substrate was examined, and it was deter- 
mined that the sputtering cross section of the adatom was enhanced with 
the increase of mass of the substrate, when all] other factors are held 
constant. This contradicts Meyerhoff's findings, but only one point was 
considered in his analyst's. 

The ratio of the sputtering yield of the bare substrate to the 
reacted substrate was evaluated for both systems. The analyses indi- 


cated that the adatoms decrease the amount of energy imparted to the 


ID 


Substrate, and reduces the sputter yield for the “hit” case. This 
momentum loss was opserved at 500 eV and e Kev. This Process could 
change at higher energy levels (greater than 2 keV) of the incident par- 
ticles. 

A significant difference was found for tungsten between the simula- 
tion and experimental and theoretical sputtering cross sections. The 
nitrogen sputtering cross section for the 500 eV bombarding energy was 
higher than for the 2 Kev energy, which was in contradiction to Winter’s 
results. A similar finding was made by Webb. This tndicates that per- 
haps the pair potentials used to model the nitrogen on tungsten are not 
correct. It has been determined that the oxygen-tungsten system behaves 
quite differently than the oxygen-molybdenum system, SO perhaps this 
trend continues for the M@PN “end W=N Systemse[ker 9 67)5 

Further simulation studies should be made to complete the compar i- 
son with Winters’ study. Studies should he made using helium and xenon 
ions in order to investigate the effect of the mass of the incident 10on. 
Additional studies should be conducted on tungsten to determine if the 
lower sputtering cross section at 2 KeV was due to nitrogen atoms being 
captured as interstitials in the tungsten lattice, or if there is a 
defect in the potential functions. Higher energy studies should be 
conducted to identify the point at which the sputtering yield of the 


substrate will be enhanced by the adatoms. 
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APPENDIX 


DERIVATION OF SPUTTERING CROSS SECTIONS 


The methodology used for the derivation of the sputtering cross 
section in this thesis 1S based on the technique used by Winters in 
meperence SS. The notation used in this appendix closely parallels that 
of Winters. 

The rate that nitrogen ts sputtered from the surface of the target 
Is given by the relation: 

CON 


-Ry = m= =2=- 0 Oy VT , (7a) 
dt 


where the following are defined: 


ime esiimrogen Sputtering rate, {atoms / cme - sec}; 
Swe= thc Mmitrogen atom comeentration at time = t, jatems / cre}: 
Oo = the sputtering cross section, fom® }; 
vt = the incident ton flux, fion / cme - Sec}, 
t = the time in seconds, {sec}; and 
ytt = the ion dose, (fluence), {ion / cme}. 


solving equation 7 yields, 


atu Onig) = -OV ty ( 8 ) 
Michie tie ttwttial mitrogen concentration at t = QO. Rearranging 
equation 8, and solving for o results in the following relation for the 


cross section; 


eM (CnN/ONg) 7 Vit. (9) 
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The Oy/Ono term 1s next examined. The following relations are 


defined, 
No = the number of nitrogen atoms on the surface at t = QO, 
NTRAJ = the number of trayectories run, also 
NTRAJ = the number of incident 10ns on the target, 
NSPUTT = the number of nitrogen atoms sputtered from the surface, 
Ao = the area of the target {cm®}. 


Consequently, for each trayectory we have 


@y = (No - NSPUTT) / Ao, 


ANd Ono = No/Ag- Thus, 


eh (N. - NSPUTT) / Ag (No. - NSPUTT) 
SNo No 7 Ao No 


Surming over all] trajectories, 





en NTRAJ (Ny - NSPUTT) , 
ONo | = 4 No 
NTRAJ 
NSPUTT | 





(No) (NTRAJ) 


NTRAJ 
1f we define »: NSPUTT - TOTSPUT, then 
a 
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On TOTSPUT 
= oo 2 ; (aie 
Ono (No) (NTRAJ) 





Incorporating the relation for the beam fluence, from equation 9 


ytelds the following relation for the cross section; 








TOTSPUTT 
oO = = in { ae Aos 

(Ny) (NTRAJ) 

TOMPSPUTT 

but = YIELD, so as a result; 
NTRAJ 
LEED 
O= - In { - Ao (eins) 


Equation 11 is the relation used to calculate the sputtering cross 
section for the adsorbed nitrogen from the tungsten and molyodenum sur - 


faces. 
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